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IN THE UNITED STATES PATENT AND TRADEMARK OFFICE ^'^^P^ ZO gg 



Attn: Special Program Law 
Office of the Office of 
the Deputy Assistant 
Commissioner for Patent 
Policy and Projects 



In re Application of: 

Ted Christopher 
Serial No: 08/746,360 
Filed: 

For: FINITE AMPLITUDE DISTORTION- 

BASED INHOMOGENEOUS PULSE 
ECHO ULTRASONIC IMAGING 

PETITION UNDER 37 C.F.R. § 1 .292 FOR THE 
INSTITUTION OF A PUBLIC USE PROCEEDING 

4 

Assistant Commissioner for PatftfjQlE.\V 
Washington, D.C. 20231 ^ 

Dear Sir: _o 

Acuson Corporation (^^Qi?y /W&^S^ petitions for the institution of a 
» public use proceeding to determine whether claims pending in the above- 

identified U.S. patent application (the "Chhstopher application") are patentable 
over the public uses described below. 

The following materials are provided in support of this petition: 
1. Two VHS videotapes labeled "Dr. Mulvahy Video II • Imaging" and 
"Mayo Clinic #CSI, Echo Date 9-28-94, Complete Copy" 

04/03/1998 BALEXAND ooooooso 08746l&°"®"'*'''®'y '^^^''^"^ *° "Mulvagh videotape"); 

^^^^'^"^ 2. ^^••W>9'hS videotape labeled "Contrast Echocardiography 

Investigations Copy A" (the "Acuson videotape"); 

3. evidentiary declarations of John M. Sheldon, Janna G. Clark, 
Paul E. Chandler, Sharon Mulvagh, M.D., Gregory L. Holley, and 
Joan C. Main; 

4. a copy of form PTO-1 449 listing the items of paragraphs 1 -3; 
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5. a check in the amount of $1 ,51 0.00. 

Additionally, Acuson hereby certifies that a duplicate copy of this petition, 
including all supporting materials, has been sent to RCT, the owner of the 
Christopher application by Federal Express on March 20, 1998, at the following 
address: 



Mr. Timothy Reckart, Esq. 
General Counsel, Secretary and 

Director of Legal Affairs 
Research Corporation Technologies 
101 North Wilmot Road, Suite 500 
Tucson, Arizona 85711-3335. 



The supporting materials listed above provide evidence of two separate 
public uses, as follows: 

The Acuson videotape shows recorded ultrasound images that were 
acquired using a fundamental transmit center frequency of 2.5 MHz (Chandler 
H 3, Main H 2) and a harmonic receive center frequency of 5.0 MHz (Chandler 
H 3, Main H 2). Selected portions of the images on the Acuson videotape show 
ultrasonic images of heart tissue obtained in the absence of contrast agent 
(Chandler HH 6-8, Holley H 7, Main Hlf 14-15), and these images were 
generated at least in measurable part as a result of second harmonic distortion 
that occurs as the transmitted ultrasound pulse propagates through the tissues 
of an animal subject. Such distortion is identified as "finite amplitude distortion" 
in Ted Christopher's January 1997 paper (Chandler H 8, Holley H 7). The 
Acuson videotape was publicly shown in the United States at the American 
College of Cardiology Scientific Sessions during the week of March 14, 1993 
(Clark H 2) and at the floor show at the American Society of Echocardiography in 
June of 1993 (Main H 12). The enclosed Acuson videotape is an accurate copy 
of the original tape (Sheldon HH 4-5, Chandler H 5 and Main H 13). 

Similarly, the Mulvagh videotape records a series of ultrasonic images 
that were acquired using a fundamental transmit center frequency of 2.5 MHz 
and a harmonic receive center frequency of 5.0 MHz (Chandler H 3, Mulvagh U 
3, Holley H 3, Main H 6). Portions of the Mulvagh videotape show resulting 



-2- 




Case No. 5050/296 

images obtained from a dog's beating heart in the absence of added contrast 
agent (Mulvagh H 4, Holley UH 4-5, Main HU 4 and 7). These portions of the 
Mulvagh videotape include images generated at least in measurable part from 
second harmonic distortion of the transmitted pulse caused by propagation 
through the tissues of the animal (Holley H 8). The unedited Mulvagh videotape 
was publicly shown in the United States on September 30, 1994 at an 
echocardiography conference entitled "Advances in Echocardiography: Contrast 
Echocardiography Perfusion Imaging" (Mulvagh HH 2 and 6, Main HH 5 and 8). 
The enclosed Mulvagh videotape is an accurate copy of the original (Mulvagh H 
5, Sheldon HH 2 and 3). 

As pointed out briefly in the preceding paragraphs, the enclosed 
declarations are evidence that ultrasound images of heart tissue, acquired at the 
second harmonic of the fundamental transmission frequency in the absence of 
added contrast agent, were in public use in the United States in March of 1993, 
June of 1993 and September of 1994. These images included a measurable 
amount of tissue harmonic signal resulting from second harmonic distortion of 
the fundamental ultrasound beam as it propagated through the tissues of the 
animal subject. Acuson wishes to have this information considered during the 
examination of the above-identified Christopher patent application, and Acuson 
submits that a public use proceeding is the proper procedure. Acuson has not 
yet reviewed the claims of the above-identified Christopher patent application, 
and is therefore not in a position to assert whether or not these claims are 
patentable over the public uses described above. Petitioner has a potential 
interest in acquiring rights in any patent maturing from the above-identified 
Christopher application, and Petitioner therefore wishes to have the above- 
described public uses considered carefully by the Patent Office in determining 
whether to grant a patent on the Christopher application. 
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Petitioner requests that the Commissioner Institute a public use 
proceeding to allow these issues to be considered. 

Respectfully submitted, 
ACUSON CORPORATION 



BRINKS HOFER GILSON & LIONE 
P.O. BOX 10395 
CHICAGO, IL 60610 
(312) 321-4200 




-William A. Webb 
Registration No. 28,277 
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Evidentiary Declaration of John M. Sheldon 



I, John M. Sheldon, declare the following facts to be correct based on my personal 
knowledge: 



1. I have been employed at Acuson since 1989 as a photo/video technician. 
Accordingly, I am very familiar with videotapes and videotape copying. 

2. I have reviewed the videotape with a thin crack in the right-hand window of ^ ^ o 
casing, referred to in Dr. Sharon Mulvagh's declaration. It is a VHS tape thartas ^ U 
excerpts that were copied from the same Betacam tape from which I made four ^ 
additional like excerpted copies on VHS tapes. Both the Betacam tape and tti^apto o 
with the thin crack include a minor copying artifact that is symptomatic of AcG§on's 
specific in-house copying set-up used when a Betacam copy is made from an original 
Super- VHS tape. I have placed my initials "JMS", and a copy number designation, 

on a label affixed to each of the four additional edited copies that I made. 

3. In addition, I made three full (unedited) copies of the Betacam tape. I have placed my 
initials "JMS", and a copy letter designation, on a label affixed to each of the three 
additional unedited copies that I made. 

4. I also made three VHS copies of an Acuson-produced S-VHS tape labeled "Contrast 
Echocardiography Investigations." I have placed my initials "JMS", and a copy 
number designation, on a label affixed to each of the three copies that I made. 

5. Finally, I made an additional four VHS copies of the Acuson-produced Betacam tape 
from which the S-VHS tape referred to in paragraph 4 was made. The Betacam tape 
is labeled "Contrast Echocardiography Investigations Loop Copy," and the label bears 
a copyright notice (with Acuson's name mis-spelled) as follows: "Copyright 1993 
Asuson". I have placed my initials "JMS" next to a copy letter designation on a 
handwritten label (entitled "Contrast Echocardiography Investigations") affixed to 
each of the four additional copies that I made. 



I hereby declare under penalty of perjury that, to the best of my knowledge, the foregoing 
is true and correct. 



Dated this 10** day of March, 1998, by 

John'T^I. Sheldon 





Evidentiary Declaration of Janna G. Clark 



CD 



I, Janna G. Clark, declare that I believe the following facts to be correct based on r^- ^ \^ 
personal knowledge: ro fj* I 

1 . I have been employed by Acuson in clinical marketing/ sales support since ^rc^ Q 
1989. O 

2. In the 1992-1993 timeframe, I worked with Paul Chandler, who was the principal 
engineer investigating contrast imaging at Acuson. In late 1992 or early 1993, we 
took a model 128 Acuson ultrasound system to Dr. DeMaria, at the University of 
Califomia at San Diego. At UCSD, under Dr. DeMaria' s direction, we imaged 
Schering's contrast agent, Levovist, in a beating animal heart. I personally did 
some of the scanning for this imaging. We recorded the images on videotape, and 
then I edited the tape for presentation. We presented it to about 100 physicians and 
sonographers on Sunday March 14, 1993, which was the eve of the American 
College of Cardiology Scientific Sessions ("ACC")- In the following several days, 
during the ACC itself, I personally presented the videotape on the show floor. I 
kept the videotape running during show hours, and answered questions posed by 
conference attendees. 

I hereby declare under penalty of peijury that the foregoing is, to the best of my 
knowledge, correct. 



Dated this 10th day of March, 1998, by 



Janna G. Clark 



Evidentiary Declaration of Paul E, Chandler 



I, Paul E. Chandler, declare that I believe the following facts to be correct based on m^personal ^ 
knowledge: 3: 

^ o 

1 . I was employed by Acuson as a Member of the Technical Staff from December.fes^o Cli 
July 1994. I was a consultant to Acuson from May 1995 to November 1997. ^gm J^ly ^ 
1994 to May 1995, 1 was neither employed by nor consulting for Acuson, due fo^ ^ O 
disability. Since approximately 1991, a significant portion of my efforts were Sftected to 
my interest in the interplay between ultrasound systems and contrast agents. 

2. In approximately 1992, at the suggestion of Reinhard Schlief, I modified an Acuson 128 
ultrasound system to support second harmonic imaging. This involved specifying the 
correct frequency pair to cause the system to transmit ultrasound energy at 2.5 MHz and 
receive ultrasound echoes at 5 MHz. To this end, I arranged to have the modification 
burned into a Programmable Array Logic (PAL) device on one of the boards, and I 
arranged to have the system software modified. I also procured a prototype wide band 
transducer from Acuson's Principal Fellow, Amin Hanafy. The whole process took a 
couple of days to complete. 

3. This first system is called "Serial Number 1 122" at Acuson. Subsequently, I instructed 
Art Tofanelli in how to modify other Acuson systems. In the next few years, pursuant to 
my requests and requests of others at Acuson, 8 additional systems were built in this 
manner. For all of these systems, when the 5.0MHz icon on the screen is highhghted, it 
means that the machine is configured to transmit ultrasound energy at a fundamental 
center frequency of 2.5MHz and receive ultrasound echoes from the anatomy and/or 
contrast agent at the second harmonic center frequency of 5MHz. 

4. I worked with Joan Main, who led the Acuson marketing contrast effort. In late 1992 or 
early 1993, Ms. Main and a cHnical marketing specialist, Janna Clark, accompanied me in 
taking Serial Number 1 122 to Dr. DeMaria, at the University of CaUfomia at San Diego. 
At UCSD, under Dr. DeMaria's direction, we imaged Schering's contrast agent, Levovist, 
in a beating animal heart. We recorded the images on videotape, and then the tape was 
edited for presentation to about 100 physicians at the ACC in 1993. I personally 
presented the edited videotape at the conference in a closed session. 

5. Today, I reviewed Acuson's Betacam tape entitled "Contrast Echocardiography 
Investigations Loop Copy," which bears a copyright nofice on the label (with Acuson's 
name mis-spelled) as follows: "Copyright 1993 Asuson". I have compared it to the S- 
VHS entitled "Contrast Echocardiography Investigations," which does not have the 
words "Loop Copy" or the copyright notice on the label. From this inspection, I am 



convinced that both tapes contain the same images and are copies of the original or copies 
of copies of the original, and that I showed one such tape at the ACC in 1993. 

6. In conducting the experiments shown on the tape, we used a protocol that included 
looking at the heart tissue at the harmonic frequency in the absence of contrast agent, and 
then looking at the heart tissue after contrast agent was added. Accordingly, the 
videotape tapes created from the experiments displayed the images at the harmonic 
frequency both before and after contrast agent was added. 

7. My purpose in running the experiments shown on the tape was to determine whether the 
second harmonic signals generated by the contrast agents could be seen above the signals 
at the second harmonic frequency known to be generated from the transmit pulse, the 
system nonlinearities, and the nonlinearities of the propagation path of the pulse through 
tissue. 



8. I have read Ted Christoper's paper, "Finite Amplitude Distortion-Based Inhomogeneous 
Pulse Echo Ultrasonic Imaging," IEEE UFFC Vol. 44(1) January 1997. The harmonic 
images on the tape that I showed in 1993 were generated at least in part as the result of 
"finite amplitude distortion" as that term is used in Ted Christopher's publication. A 
copy of this publication is attached hereto. 
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Finite Amplitude Distortion-Based 
Inhomogeneous 
Pulse Echo Ultrasonic Imaging 



Ted Christopher 



Abstract—Ultrasonic pulse echo imaging in inhomoge- 
neous media suffers from significant lateral and contrast 
resolution losses due to the defocusing effects of the inho- 
mogeneities. The losses in lateral and contrast resolution 
are associated with increases in the width of the main- 
beam and increases in sidelobe levels, respectively. These 
two forms of resolution loss represent significant hurdles to 
improving the clinical utility of biomedical ultrasonic imag- 
ing. A number of research efforts are currently under way 
to investigate the defocusing effects of tissue and to con- 
sider corrective measures- All of these efforts assume \\n- 
ear propagation, and base the image-formation process on 
the reception of the transmitted pulse. A novel pulse echo 
imaging scheme in which the image is formed using the fi- 
nite amplitude distortion components of the received pulse 
is considered here. Alternatively, this could be described 
as image formation using the nonlinearly-generated higher 
harmonics. 

In homogeneous beam propagations, it has been estab- 
lished that the sidelobes of nonlinearly-generated higher 
harmonics are much lower than their linear counterparts. 
Computations considered here suggest that this relation- 
ship also holds for the case of propagations through ab- 
dominal wall and breast wall tissue. These computations 
also suggest that the lateral resolution limits imposed by 
abdominal wall and breast wall tissue are slightly smaller 
for nonlinearly-generated second harmonics than for their 
linear counterparts. The resulting potential of these higher 
harmonics to improve image resolution is investigated. 



I. Introduction 

THE LOSS OF CONTRAST and lateral resolution due to 
phase aberrations produced by inhomogeneities in 
tissue is a significant problem in biomedical ultrasonic 
imaging [ll-[l3]. A variety of pliasc aberration cor- 
rection schemes have been considered [2]45]-[7],[8],[10l- 
[1114131. These schemes have thus far demonstrated lim- 
ited clinically-relevant success. 

A new approach, and one independent of the existing 
approaches, involves using the fmilc amplitude or nonlin- 
ear distortion of the propagating beam to form an image. 
In this approach, the ultrtusound image would be based 
on the received finite amplitude distortion components (or 
nonhnearly-generated higher harmonics) associated with 
the transmitted signal. In the simplest case in which the 

M.-xnuscripi received March 5, 1996: accepic-d -luly 22, 199G. 

T Christopher is with the DcparinicuL of r::icclrical l£n- 
ginocring aiul Rochester Cenler for Biomedical Ultra.sound 
University of Rochester. Rochester, New York 14627 (e-inail: 
chrislo'^i'oe.rochesier.echi). 



transducer emits negligible energy in the second harmonic 
bandwidth, such an image could he formed by adding an 
initial high pass filtering of the received signal. In gen- 
eral, such an image could be formed by using a two pulse 
transmit, receive, normalize, and then high pass filtering 
scheme. Such a two pulse schenic could be used to remove 
the source or linear content in the second and higher har- 
monic band widths. 

In the following section, the finite amplitude produc- 
tion of higher harmonic beams will be considered for the 
case of an apodized and focused Gaussian transducer op- 
erating in a non-aberrating liver medium. Both contmu- 
ous and pulsed fields will be considered. The associated 
imaging-relevant sidelobe levels and mainbeam widths of 
the second and third harmonics will be considered there. 
In Section III, the requisite field amplitudes and their re- 
lationship to the proposed mechanical index (MI) will be 
considered. 

In Section IV, the essential nuithematics describing the 
formation of higher harmonic beams will be presented. 
This presentation provides insight into the homogeneous 
propagation results shown in Section II and a basis for un- 
derstanding the subsequent inhomogeneous propagations 
considered. In Section V, inhomogeneities are added to 
the previously considered liver path propagation example. 
The inhomogeneities are in the form of planes of time de- 
lays, which were measured for abdominal wall and breast 
wall layers [4l,[7]. Of particular interest in this section are 
the results depicting the relative impacts to the mainlobe 
widths and sidelobe levels of linear (fundamental) beams 
and nonlinear (higher harmonic) beams. Finally, in Sec- 
tion VI. the prospects for nonlinear imaging are discussed. 



II. Homogeneous Nonlinear Propagation 
Examples and Analysis 

A. Continuous Wave Soiltccs 

Several computational models exist which can accu- 
rately describe the finite amphtude propagation of a con- 
tinuous beam. Such models can be extended to compute 
the case of a propagating pulse, ^us well. These models ac- 
count for the effects of diffraction, nonlinearity, and ab- 
sorption. The excellent accuracy possible with such mod- 
els is depicted in [U]-\\7l The NLP model [16] was used 
to compute the following linear and nonlinear examples. 
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Fig. 2. The harmonic velocities used in propagating the 2 MHz Gaus- 
sian beam through hver. 
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Fig. 1. Linear propagation results for a focused 2 MHz Gaussian 
transducer operating in a Mver medium, (a) The source plane ampli- 
tude profile, (b) The on-axis amplitude, (c) The radial focal plane (z 
= 6 cm) beam profile. 



This model has been updated to account for the effects of 
dispersion. The associated harmonic velocities were com- 
puted using a novel algorithm [18]. The NLP model also 
can account for the effects of refraction and reflection at 
a planar fluid boundary which, though not utilized in the 
following examples, could be relevant to modeling some in 
vivo propagations. 

All of the following propagations are for a Gaussian 
apodized, axially symmetric focused source. This form of 
device should offer excellent image quality and produces a 
field comparable to the two-dimensional array-based trans- 
ducers now being developed. The phase aberrated fields 
of these devices, considered in Section V, were not ax- 
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Fig, 3. Nonlinear propagation results for the focused 2 MHz Gaussian 
transducer. The propagation medium here and in all subsequent fig- 
ures is liver, (a) The on-axis amplitudes of the fundamental, second, 
and third harmonics (top to bottom curves, respectively), (b) The 
focal (z = (y cm) linear (solid curve) and nonlinear (dashed curve) 
waveforms. 
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FiK 4 Fundamental and second harmonic, one- and twc^way focal 
rig. t. ruuu» craled norma zed one-way focal 

plane beam patterns, (a) The log-scaieQ, norm 
Plane profiles of the 2 MHz fundamental (soUd curve) the 4 MHz 
sto d' harmonic (long dash curve), and the 4 MHz fundamental 
(short dash curve), (b) The corresponding tw<^way profiles. 

ially symmetric. In order to model these fields a two- 
dimensional version of NLP was utilized. 

Many biomedical imaging devices are not axially sym- 
metric though. The finite amplitude beams produced by 
such devices should be well described by the computations 
for the axially symmetric transducer considered here. In 
measuring the nonlinear harmonic generation rom an un- 
focused rectangular source. Kamakura, Tarn, Kumamoto, 
and Ueda [19) noted that "the [nonlinear] harmonic pres- 
sure levels in the far field (were) almost the same as fronri 
a circular aperture source with equal f^^^^^^^f 
initial pressure, independent of the source levels. Though 
this result was obtained for only one device with a ra- 
tio of source side lengths or aspect ratio of 11 to 6 it 
should be the, case that the higher harmonic pressure lev- 
els associated with a non-axially symmetric device should 
be approximately equal to those of the corresponding ax- 
ially symmetric source. More importantly, the lack of ax- 
ial symmetry should not affect the relative side obe ad- 
vantages exhibited by the nonlinearly-generated harmon- 
ics in a homogeneous propagation. Experimental results 
[201 support this point. These homogeneous path sidelobe 
level advantages are the b^isis for the imaging-relevant ad- 
vantages in an inhomogeneous propagation. Ihe lesults 
displayed in Fig. 1 were obtained by <^°'^P"t'"S MH^ 
ear, liver path propagation of the field of a focused 2 MHz 
Gaussian source using the NLP beam propagation model. 




• 

1 5 2 2.5 

TIME (microsec) 

Fig 5 (a) The on-axis source plane (dashed curve) and suhs<.^^en, 
focal plane (solid curve) 2 MHz-centered pulses. The per.od of the 
pulses' s 8 microseconds, (b) The corresponding SP-^-'--,^ J^*; 
focal nonlinear distortion pulse obtained by =°"f V^^^" 
form using only the focal spectral information shown in (b) from 3 

to 8 MHz. 



The NLP model propagates a planar, normal velocity de- 
scription of the acoustic field. No inhomogeneit.es or phase 
aberrations were accounted for in this propagation or any 
of the subsequent propagations considered in this sect.oii. 
The relevant liver propagation parameters used were c - 
1570 m/s, p = 1.05 g/cm^ a = 0.03 Np/cm, and b = L3 
(where a and 6 are the coefiicients describing absorption 
in a power law form) (21), (22). 

The Gaussian amplitude shading of the source plane 
normal velocity field w.vs such that the half-ampl. ude^j. 
dial distance was 0.84 cm. The on-axis, source plane RMS 
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Kig. 6. Focal plane profiles from ilie 2 MHz continuous wave propaga- 
tion (of Section II. .A.) and the 2 MHz-centered pulse propagation (of 
Section Il.B.). The sohd curve depicts the 4 MHz second harmonic 
profile. The long dashed curve depicts the temporal peak amplitude 
profile computed using the frequency bandwidth of 3 to 8 MHz. The 
short dashed curve depicts the corresponding peak amplitude profile 
for the frequency range of 3 to 5 MHz. 



acoustic intensity (pc|up/2, where u is the acoustic par- 
ticle velocity) for this field was 2 VV/cm^. The radial ex- 
tent of the source was 1.5 cm. The source plane field was 
focused using a spherically-focusing factor {e^^^''\ where 
e{r) = (27r//c)\/r2 + F^). The geometric focal length F 
was 6 cm and the sound speed (c) used to compute 9{r) 
was that of water (1500 m/s). 

Figs, la and (b) depict the normal velocity magni- 
tudes of the Gaussian transducer's source plane and on- 
axis fields, respectively. Fig. Ic displays the focal plane {z 
~ 6 cm) radial profile of the 2 MHz field. The drop in the 
magnitude of the field from the mainlobe to the first side- 
lobe in Fig. Ic is 36 dB. In the absence of strong medium 
phase aberration, this should allow the device to produce 
high contrast images. 

The same 2 MHz Gaussian source was then propagated 
nonlinearly through the same liver path. The nonlinear pa- 
rameter p used to represent liver was 4.7 [23]. The NLP 
model used 4 harmonics (2. 4. 6. and 8 MHz) to com- 
pute the pre- focal region (~ = 0 to z = 3 cm) propagation 
and up to 10 harmonics to represent the subsequent focal 
and post- focal region propagation. The harmonic veloci- 
ties were computed using a novel algorithm described in 
[18]. The fundamental or 2 MHz component had a propa- 
gation speed of 0.157 cm/microsecond (given above as c). 
The discrete harmonic velocities used by NLP to compute 
the 2 MHz propagation are shown in Fig. 2. 

Fig. 3a displays the axial magnitudes of the fundamen- 
tal, second harmonic, and diiixl harmonic fields as com- 
puted for the nonlinear propagation. The fundamental or 
2 MHz axial curve is only slightly different from the corre- 
sponding linear curve shown in Fig. lb. At 2 = 6 cm, the 
lionlinear 2 MHz curve is about 1% lower than the 2 MHz 
linear curve. This difference was due to growth of the 
higher harmonics in the nonlinear propagation. In Fig. 3b 
the corresponding focal (z = 6 cni) pressure waveforms 
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Fig. 7. (a) Log-scaled first, second, and third harmonic axial ampli- 
tudes (top to bottom, respectively) for the focused 2 MHz Gaussian 
transducer, (b) The corresponding log-scaled, focal plane radial beam 
profiles. 
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Fi- 9 Homogenous and abdo.ninal wall-jittered tissue path beam propagation results for the 2 and 4 MHz focused Gaussian transducers. 
Focal plane (one-way) profiles, the solid curves are the un-jittered or homogeneous tissue path results, (a) 2 MHz fundamental diameters, 
(b) Average radii from the 2 MHz focal plane profiles, (c) 4 MHz second harmonic diameters, (d) Average radii from the second harmonic 
-I MHz profiles, (e) and (f) same format for the 4 MHz fundamental profiles. 



are displayed. The pressure wa\'eforms were obtained by 
converting NLP's normal velocity output to pressure using 
the impedance relation (all subsequent pressure waveforms 
were obtained in this way). 

The 2 MHz or fundamental beam pattern associated 
with the nonlinear propagation is almost identical with the 
2 MHz beam pattern of the linear propagation. Only in a 
linearly scaled overlay plot of the two beam patterns are 
there visible differences. These differences are very small 
and are limited to the near axis portions of the beam pat- 
terns. Only at much higher source intensities [16], [24) are 
the effects of nonlinearity significant to the details oi the 
fundamentars field. These results are consistent with the 
empirical observation tliat linear modeling of biomedical 



ultrasonic device fields accurately describes their (linear- 
based, homogeneous path) imaging performance. 

Fig. 4a depicts the 2 MHz fundamental and 4 MHz sec- 
ond harmonic focal plane amplitude profiles. Also shown 
in Fig. 4a is the corresponding 4 MHz fundamental profile. 
The 4 MHz fundamental result was obtained by comput- 
ing the linear propagation of the same Gaussian trans- 
ducer operating at a source frequency of 4 MHz. All three 
beam profiles have been normalized to have on-axis field 
amplitudes of 1. The finite amplitude distortion-generated 
second harmonic focal profile has a slightly broader main- 
lobe than the corresponding (4 MHz) fundamental profile. 
The radial half-ami)litu(le distance of the second harmonic 
profile is 36% greater tlian tliat of the 4 MHz fundamental 
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Fig. 10. Additional results from the abdominal wall-jittered propagation path case shown in Fig. 9. (a) 2 MHz fundamental, 4 MHz second 
harmonic, and 6 MHz third harmonic average two-way radii, (b) 2 MHz, 4 MHz, and 6 MHz fundamental or linear average two-way radii, 
(c) 4 MHz second harmonic and fundamental average two-way radii, (d) 6 MHz third harmonic and fundamental average two radii. The 4 
and 6 MHz fundamental results in (c) and (d) are the shorter dash curves. 



profile (0.0983 cm versus 0.0723 cm). The second harmonic 
profile also has much lower sidelobes than the 4 MHz fun- 
damental profile. 

For imaging purposes, the two-way focal plane beam 
pattern of the Gaussian transducer is of interest. The two- 
way focal beam pattern should account for both the char- 
acteristics of the transmitted pulse in the focal plane and 
the corresponding sensitivity of the transducer to pulses 
reflected back from this plane. For linear propagations, 
the two-way beam pattern for a given depth can be ob- 
tained by squaring the corresponding transmit or one-way 
beam pattern. Fig. 4b depicts the normalized two-way lin- 
ear beam patterns for the Gaussian transducer operating 
at 2 and 4 MHz. These curves were obtained by squar- 
ing the corresponding one-way or transmit beam patterns 
shown in Fig. 4a. 

Also shown in Fig. 4b is the two-way focal plane beam 
pattern associated with the 4 MHz second harmonic field. 
Since the amplitudes of the reflected pulses are much 
smaller than the transmitted pulses, the propagation of the 
reflected field back to the transducer is essentially linear. 
Thus the two-way, second harmonic focal plane beam pat- 
tern was obtained by multiplying the corresponding one- 
way pattern shown in Fig. 4a by the 4 MHz fundamental 
one-way pattern also shown in Fig. 4a- 

The second harmonic's two-way beam pattern has a 
half-amplitude mainlobe width (or -6 dB beamwidth) 
that is 12% greater than that of the corresponding fun- 



damental beam pattern. The -20 dB beamwidth of the 
second harmonic is 13% greater than that of the 4 MHz 
fundamental. The sidelobe advantage displayed in the fo- 
cal plane profiles of Fig. 4a is maintained in the two-way 
results of Fig. 4b. 

These homogeneous liver path results suggest that the 
second harmonic field of a focused, apodized transducer 
might oflter advantages in contrast resolution over the cor- 
responding or same-frequency transducer field. The results 
also suggest that these advantages would be slightly offset 
by losses in lateral resolution. 

B- Pulse Sources 

A pulse propagation was next considered for the 2 MHz 
Gaussian source. The on- axis, source plane pressure pulse 
used is displayed in Fig. 5a as the dashed curve. This pulse 
was computed by applying a Gaussian window to a 2 MHz 
cosine. The peak pressure of the pulse was the same as for 
the 2 W/cm"^ continuous case considered in Section II. A. 
The apodization and the spherical focusing of the source 
field were also the same as in the previous continuous wave 
case. The initial source plane pulse consisted of 128 sam- 
ples across 8 microseconds. 

The magnitude of the Fourier transform of the 8- 
microsecond source pulse is shown in Fig. 5b. A straight- 
forward implementation of the nonlinear imaging idea re- 
quires that there be negligible overlap between the source's 
spectral bandwidth and that of the nonlinear second har- 
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Fig. 11. Two-way averaged radial results from five abdominal wall- 
jittered propagations, (a) Normalized two-way average radii for the 
2 MHz fundamentals (solid curve), 4 MHz second harmonics (long 
dash curve), and 4 MHz fundamentals (short dash curve), (b) The 
corresponding radially integrated magnitudes. These curves have 
been normalized to have a unity value at r = 1.2 cm. 



monic. (More generally, this would also insure negligible 
overlap between any of the successive harmonic spectral 
bands.) The source spectrum shown in Fig. 5b meets this 
requirement. Nonlinear images based on source pulses with 
broader spectrums or with significantly more energy in 
the second harmonic bandwidth than the one depicted in 
Fig. 5 could be obtained by using an alternative nonlinear 
imaging scheme described in Section VI. 

The source plane was then defined using the 64 har- 
monic Fourier transform of the source pulses. This multi- 
harmonic source radius was then input into a pulse- 
propagating version of the NLP model (a model similar 
to the lithotripter model [25]). The focal output of the re- 
sulting nonlinear pulse propagation is also shown in Fig. 5a 
and 5b as solid curves. The focal pulse waveform has 
slightly smaller peak positive and negative pressures than 
the corresponding continuous waveform shown in Fig. 3b. 
Consistent with the smaller amplitudes, the focal pulse is 
also less distorted than the continuous waveform. The ra- 
tio of the second harmonic's focal amplitude to that of the 
fundamental's for this pulse propagation was 70% of the 
same ratio for the corresponding continuous source con- 
sidered in Section ILA. 

Fig. 5c shows the waveform associated with only the spec- 
tral bands of the second, third, and the first half of the 
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Fig. 12. Two-way averaged radial results from five abdominal wail- 
jittered propagations, (a) Normalized two-way average raxiii for the 
4 MHz fundamentals (solid curve), 8 MHz second harmonics (long 
dash curve), and 8 MHz fundamentals (short dash curve), (b) The 
corresponding radially integrated, normalized magnitudes. 



fourth harmonic (3 through 8 MHz), This waveform was 
computed as a high pass filtered reconstruction of the spec- 
tral information depicted in Fig. 5b. A rectangular window 
with a transition at 3 MHz was used in filtering the trans- 
form data. The straightforward nonlinear imaging approach 
considered here would use a distortion pulse like that shown 
in Fig. 5c in order to image the scattering medium. 

Not shown for the pulse propagation considered are the 
axial and radial harmonic descriptions. These descriptions 
were found to be identical in form to those computed for 
the nonlinear propagation of the continuous 2 MHz field. 
All other Gaussian-windowed cosine pulses were found to 
produce identical axial and radial harmonic patterns. 

The temporal peak amplitude profile of the high pass 
filtered- focal plane data also very closely followed the form 
of the 4 MHz second harmonic's amplitude profile. The 
temporal peak amplitude profile for the pulse propagation 
considered in this section is shown in Fig. 6. This profile 
wtAS computed using the focal plane frequency data in the 
range of 3 to 8 MHz (like the waveform in Fig. 5c). Also 
shown in Fig. 6 are the corresponding 4 MHz second har- 
monic profile and a second harmonic bandpass filtered (3 
to 5 MHz)-temporal peak amplitude profile for the same 
pulse propagation. The similarity of the two peak ampli- 
tude profiles suggests that there is very little energy above 
5 MHz for this pulse propagation. 
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Fig. 13. Two-way averaged radial results from five breast-jittered 
propagations, (a) Normalized two-way average radii for the 2 MHz 
fundamentals (solid curve). 4 MHz second harmonics (long dash 
curve), and 4 MHz fundamentals (short dash curve), (b) The cor- 
responding radially integrated magnitudes. These curves have been 
normalized to have a unity value at r = 1.2 cm. 



As shown for the second harmonic bandwidth in this 
case, it appears that the higher harmonic peak amphtude 
patterns produced by any reasonable source pulse can be 
computed using an appropriate continuous approximation 
of the source. This is a very convenient relationship. 



Ill- Requisite Field Amplitudes 

In order for the nonlinearly-generated higher harmonics 
to be available for imaging in an inhomogeneous media, the 
received higher harmonics amplitudes can't be too low rel- 
ative to the received linear (or transmitted) signal and the 
receiver's dynamic range. Also, to be available for imag- 
ing use, the associated in vivo field amplitudes must not 
correspond to a mechanical index greater than 1.9. 

For the previously considered 2 MHz liver-path prop- 
agation, the amplitudes at the focus {z - 6 cm) of the 
fundamental and second harmonic were 0.943 MPa and 
0.166 MPa, respectively. The second harmonic amplitude 
was thus 15.09 dB below the amplitude of the 2 MHz trans- 
mitted wave. As a wave reflected at the focal point travels 
the 6 cm back to the transducer/receiver, frequency de- 
pendent attenuation reduces the second harmonic by an 
additional 5.62 dB relative to the 2 MHz component. The 



received 4 MHz component would thus bo 20.71 dB bck)w 
the transmitted 2 MHz component, (ignored here is tlic 
possible additional relative weakening of the second har- 
monic received .signal due to coherent reflectors at the fo- 
cus in combination with the smaller mainlobe of the second 
harmonic. This effect should not l^e significant to many 
bio-ultrasound imaging applications, though.) The corre- 
sponding figures for tlie 6 MHz third harmonic were 27 dB 
down at the focus and 39.23 dB clown for received signals. 

Table I displays the second and third harmonic received 
levels corresponding to on-axis, source plane intensity val- 
ues of 0.5, 1, 2, 4, and 8 W/cm^, For a given source plane 
intensity, pulsed devices would have slightly larger received 
level differentials than those shown in Table 1. For the 
pulsed propagation considered in the previous section, this 
additional gap would be 0.86 dB for the second harmonic 
(based on a second harmonic bandpass reconstruction and 
a comparison of the received peak positive pressure levels). 

Current biomedical ultrasonic imaging transducers have 
dynamic ranges of about 100 dB. Even with decreased sen- 
sitivity above the transmit frequencies, these devices may 
be capable of creating second harmonic images. This ca- 
pability may well have been demonstrated by the creation 
of second harmonic contrast agent-response images. Alter- 
natively, a separate receiver device with appropriate fre- 
quency response in the desired nonlinear distortion band- 
width could be used. 

The effect of focal length on the received second har- 
monic levels for this Gaussian transducer operating at 2 
W/cm^ is shown in Table IL From a focal length of 4 cm 
up to a focal length of 12 cm, the received second harmonic 
levels dropped off by 7.32 dB. Also shown in Table II are 
the corresponding focal second harmonic levels. These lev- 
els remained very constant and thus revealed the decreases 
in received levels as almost entirely due to increased return 
trip distances- 

The peak positive and negative pressures of the pre- 
dicted in vivo nonlinear waveform shown in Fig. 3b were 
1.12 and -0.84 MPa, respectively. The -0.84 peak rarefac- 
tion pressure corresponds to a mechanical index (MI) of 
0.59. The highest acceptable level for the mechanical index 
is 1.9. Shown in Table I are values of the minimum focal 
pressure and associated mechanical index for this and four 
other values of source plane intensity. Note that the min- 
imum pressures and thus mechanical indices given in Ta- 
ble I have been corrected for the effects of noniinearity. A 
linear-only computation would have predicted larger neg- 
ative pressures and MI values, in particular at the highest 
two source intensity levels. 

The numbers shown in Table I suggest that for in vivo 
propagations similar to the one considered here, finite am- 
plitude distortion-based images should be readily obtain- 
able within the current mechanical index safety limit. Even 
at the lowest source intensity case considered (0.5 VV/cm"). 
a largely second harmonic-based image could be obtained 
by simply filtering out the transmitted freque!icy(ies). Ad- 
ditionally, the second and third harmonic received levels 
could offer some real-time feedback on the magnitude oi 
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TABLI^I 

ON-AXIS (SPAT.AL PKAK. SOI HCK PLANK 1NTKNS1TV M-RSlS I.KC K I VKD SKCOND AND 1 liim) HARMONIC LKX KI.S 
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On-axis source intensity Received 2nd Received 3rd Mininn.m pressure at the j^;^ ^^ 

(RMS w/cn.-M Imrnionic level (dB) liarnionic level (dB) geometric focus (MPa) ('""''l/l/V j 

5^;^^^ 50^8 -0^^ 

l' 23.64 45.02 -0.61 0.43 

.y 90 73 39.23 -0.84 0.59 

17 94 33.66 -1 13 0.80 

8 15.37 28.59 -1^0 \m 
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the focal field amplituderi themselves. Finally, the results 
di.si)layed in Table II suggest that second harmonic imag- 
ing may be available at a wide range of focal depths. 



IV. Basics of Nonlinear Beam Formation 

The formation of the higher harmonic constituent 
beams in a propagating finite amplitude beam is a contnv 
uous process. In the case of the 2 MHz Gaussian-shaded, 
focused beam considered thus far, the 4 MHz second har- 
monic, the 6 MHz third harmonic, and additional higher 
harmonics are continuously and cumulatively produced by 
the beam as it propagates away from the source. Of in- 
terest here is the production and focusing of these nonlin- 
ear higher harmonic beams between the source and focal 
plane. 

The origin of the higher harmonic beams is the ongoing 
nonlinear distortion of the propagating waves comprising 
the (total harmonic) focused beam. The physical effects of 
diffraction and absorption concurrently act on the higher 
harmonic beams and thus further define their propagation 
as well as contribute to changes in the resulting focused 
beatn. The NLP model assumes that the nonlinear or finite 
amplitude distortion acts in a plane wave fashion on the 
waves comprising the focused beam [16). The NLP model 
uses the frequency domain solution to Burgers' equation in 
an incremental fashion to account for this plane wave 
distortion approximation. (See [IG] for further details in- 
cluding the inclusion of the local directionality of the beam 
front, which is neglected in the following presentation.) 

The frequency domain solution to Burgers' equation as 
used in the NLP model can be written 

n= 1.2.... ..V ^ (1) 

where / is the fundamental frequency and Un is the /?th 
term in an N term complex Fourier series describing the 
temporal normal velocity waveform at a point in the radial 
description of the field. The first summation term in the 
parentheses represents the accretion of the nth harmomc 
by nonlinear combination of other harmonics that have a 



TABLE II 

Focal and tuxkived sfxono uahmonic levels (reu-viivk to the 
fcnda.\ie.\tal) versus focal length. The geometric foc:ls of 

THE 2 MHz GaUSSL\N TRANSDICER AND THE PRO PAG. ATI ON LENGTH 
of the model were set to THE FIVE DISTANCES (6 CM IN THE 
PREVIOL'SLV CONSIDERED CASE). THE ON-AXIS SOCRCE INTENSITY 
WAS 2 W/CM^ FOR ALL OF THE PROPAGATIONS. 



Focal length 
(cm) 


2nd harmonic level 
at the focus (dB) 


Received 2nd 
harmonic level (dB) 


4 


15.99 


19.76 


6 


15-11 


20.73 


8 


15.01 


22.53 


10 


15.29 


24.69 


12 


15.83 


27.07 



sum frequency of nf. The second summation term (with 
conjugation) represents the depletion of the nth harmonic 
to other harmonics with a difference frequency of nf. For 
the case of a Fourier representation of a (periodic) pulse 
waveform, this accretion and depletion of harmonics re- 
sults in some interesting phenomenon, including the pro- 
duction of a distortion bandwidth below the fundamental 
bandwidth. 

Of interest to the results presented in this paper are the 
terms in (I) contributing to the growth of the second har- 
monic and, to a lesser extent, the third harmonic. When 
n equals 2 in (1) the positive contributions to ^ come 
from the Iil\Ui term comprising the first sumnuition. The 
negative contributions to the second harmonic are repre- 
sented by the second summation in (1) and, for the prop- 
agations under consideration here, can be approximated 
by the second term in that summation, 2tx3iiJ. Even this 
term, though, is negligible throughout most of the source 
to focal region propagation due to the relatively small am- 
plitudes of the third harmonic. Thus, the nonlinear pro- 
duction of the second harmonic throughout most of the 
relevant propagation region is simply proportional to the 
square of the amplitude of the fimdamental. 

When n equals 3 in (1), the positive contributions to 
^ come from tlie first two terms comprising the first 
sununation. These two terms sum to 3uiU2- The first term 
of the corresponding negative contributions to the thud 
harmonic is 'hnu\. This term and its successors are ueg- 
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ligible for all but the focal region of the highest intensity 
propagations considered here. 

The harmonic sources of the second harmonic and third 
harmonic are thus \uiUi and 3uiU2, respectively. At a 
given point in the field of the propagating beam then, 
the finite amplitude production of the second harmonic is 
proportional to the square of the fundamental harmonic's 
amplitude. The production of the second harmonic off the 
beam axis is very small, since the ampHtude of the funda- 
mental beam there is quite small. The third harmonic is 
produced in proportion to the product of the first and sec- 
ond harmonics, and thus its nonlinear production is even 
more strongly weighted towards the beam axis. 

Also of interest is how the production of the second and 
third harmonics varies with the z coordinate. Neglecting 
the effects of absorption and approximating the effects of 
focusing by assuming spherically-converging wave propa- 
gation, we should expect the amplitudes of the fundamen- 
tal harmonic's mainbeam to increase linearly with distance 
from the source to the focus. At z = F/2, the on-axis 
ampHtudes of the beam should be about twice the cor- 
responding source plane ampHtudes. Following from the 
same approximations, the fundamental's mainbeam width 
at 2 = F/2 should be about half its source plane width. 
Thus at 2 - F/2, the area of the fundamental's mainbeam 
should be about 1/4 the corresponding source plane area. 
Together this relation and the previous amplitude relation 
suggest that the second harmonic beam production rate 
versus z is constant {z < F), with rate losses due to the 
diminishing fundamental mainbeam area balanced by the 
concurrent gains due to increased fundamental amplitudes. 
In spite of the gross approximations used in deriving this 
relation, initial computational analysis suggests that it is 
a useful approximation. 

An analogous derivation for the relationship between z 
distance and the nonlinear production of the third har- 
monic will not be pursued here. The fact that the third 
harmonic production is proportional to the product of the 
fundamental and second harmonic amplitude, though, im- 
plies that the production of the third harmonic is strongly 
weighted towards the focal region. Fig. 7a displays the 
log-scaled, axial velocity amplitudes of the fundamental, 
second, and third harmonics. The source was the same fo- 
cused 2 MHz Gaussian source considered in Section ILA. 
The medium's parameters were again those of liver. Con- 
sistent with the above observations, the amplitudes of the 
second harmonic exhibit a relatively large "gain" in its 
growth from low near field values to significant focal am- 
plitudes. The third harmonic exhibits an even higher gain, 
roughly duplicating the growth in gain from the fundamen- 
tal to second harmonic. Both harmonics, though, display 
post-focal region amplitude declines which parallel those 
of the fundamental. 

Fig. 7b displays the corresponding log-scaled, focal 
plane (2 = 6 cm) radial beam profiles. The fundamental 
profile drops 49. S dB over the 1 cm radial range shown. The 
second harmonic roughly squares this decline in dropping 
85.3 dB. The third harmonic then continues the relation- 



ship in dropping 120.5 dB. These declines reflect the sec- 
ond and third harmonic, finite amplitude production rates 
discussed above. At 2 = 8 cm this relationship between 
the harmonic beam profiles continued to hold. 

In Fig. 8 the axial amplitudes of the 4 MHz second har- 
monic beam are overlaid with those of the corresponding 
4 MHz fundamental harmonic beam. The on-axis source in- 
tensity of both the respective propagations was 2 VV/cm^. 
The axial curves have been normalized to be unity at 2 — 6 
cm and log-scaled. The 4 MHz second harmonic focal ampli- 
tude was originally (pre-normalization) 15.0 dB below the 
corresponding 4 MHz fundamental value. The two curves 
are close through the focal region, then depart shortly after 
2 — 6 cm as the 4 MHz linear curve rapidly declines. 

The results displayed in Figs. 7 and 8 suggest that the 
4 MHz second harmonic beam could be less susceptible to 
the defocusing effects of near field phase aberrations than 
a 4 MHz fundamental beam. Since only a fraction of the 
second harmonic beam forms in the near-field, only this 
fraction could be redirected or defocused by near-field jit- 
ter. The corresponding 2 MHz fundamental beam, though, 
would pass in its entirety through the aberration and suf- 
fer the consequent defocusing effects, including increased 
focal plane sidelobe levels. Secondarily, these higher 2 MHz 
sidelobe levels could in turn increase the off-euxis nonlinear 
production of the 4 MHz second harmonic. The next sec- 
tion will consider the effects of phase aberrations on linear 
and nonlinear harmonic beams. The potential of the non- 
linear harmonics to improve on the inhomogeneous imag- 
ing performance of linear harmonics will depend on their 
ability to maintain lower sidelobe levels. 

V. The Effects of Phase Aberrating 
Inhomogeneities on Linear and 
Nonlinear Beam Propagations 

In order to investigate the effects of tissue-based phase 
aberration on the characteristics of linear and nonlinear 
beams, planes of phase delay or jitter were introduced into 
linear and nonlinear propagations of the focused Gaussian 
transducer. These phase delay planes were computed us- 
ing measured time delays from five abdominal wall layers 
[4) and five breast wall layers [7]. The five abdominal wall 
specimens had layer thicknesses of 2.5, 2.0, 1.5, 1.5, and 
1.0-3.0 cm (a non-uniform layer), giving an average thick- 
ness of 1.9 cm. The five breast wall specimens had layer 
thicknesses of 1.5-2.5, 3.0-3.5, 3.5, 4.0, and 2.0-2.5 cm, 
giving an average thickness of 3.0 cm. The average thick- 
nesses of the non- uniform layers were used to compute the 
five-layer averages. All of the abdominal wall and breast 
wall layers contained an outer skin layer. 

Each of the mciisured abdominal wall time delay planes 
was converted to an equivalent 2 MHz phase delay plane. 
Each of these delay planes was then scaled by 0.5 (i.e., 
each phase delay w^\i; reduced by a factor of 2) and then 
applied twice to a given beam propagation to represent 
the cumulative aspect of the actual tissue delays. In ap- 
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plying n single delay plane, the 2 MHz phase delay values 
were scriled for appropriate apj^licatiou to each harmonic 
present in the computed field. The first delay plane was 
encountered by the propagating field at z = 0.5 cm and 
the second plane at 2 = 1.5 cm. Further, subdivision (and 
subsequent application) of the abdominal wall delay data 
did not appear to be necessary since it did not significantly 
change the resulting focal plane fields. Thus, it appeared 
that the 2 delay plane application scheme satisfied the thin 
lens api^roximation [26). The breast delay planes were like- 
wise applied in two steps, the first at 2 ^ 1 cm and the 
second at z ^ 2 cm. For both tissue types, the z place- 
ment was selected to represent the average slice thickness 
and also be convenient for the A: step size utilized by the 
linear propagation. 

Fig. 9 (a-f) depicts pairs of one-way focal plane har- 
monic amplitude diameters and the corresponding aver- 
age radii for an unjittered (or homogeneous) path and an 
abdominal wall-jittered propagation path. The propaga- 
tion parameters of both media were again those of liver. 
In Fig. 9a the corresponding focal plane diameters for the 
unjittered and jittered 2 MHz fields are overlaid. In Fig. 9b 
the corresponding average radii are shown for the 2 MHz 
fields. The average radii were obtained by averaging the 
focal plane grid of amplitudes around the axis. The corre- 
sponding results for the second harmonic 4 MHz and the 
fundamental 4 MHz fields are shown in Figs. 9(c-d) and 
9(e-f): respectively. 

In Fig. 9, two basic eflPects of jitter are visible. The first is 
the increased sidelobe levels associated with the defocusing 
of the abdominal wall phase delays. This sidelobe effect 
is more prominent for both of the 4 MHz fields than for 
the 2 MHz fundamental field. The shorter wavelengths of 
these 4 MHz fields allows for greater de-focusing by the 
phase screen. The second harmonic 4 MHz field also gets 
a sidelobe level increase from the corresponding increase 
in the 2 MHz field's sidelobe levels. 

The second impact revealed in Fig. 9 is that there are 
small changes in the mainlobes of all three harmonic fields. 
In particular, in the diameter figures a shift in the peak or 
center of the jittered mainlobes can be seen. In the corre- 
sponding average figures, though, the impact of the jitter 
appears to be negligible, down to approximately 20 dB 
below the peak on-axis value. Thus, the impact of this ab- 
dominal wall-jittering doesn^t appear to involve any sig- 
nificant broadening of these mainlobes. Finally, both of 
the jittered 4 MHz mainlobes show decreases in peak am- 
plitudes which are consistent with the increased energy 
present in the sidelobe regions. 

Fig. 10 shows two-way average radii results for tlie same 
abdominal wall-jittered propagations. Each of the curves 
shown was obtained by radially averaging the correspond- 
ing two-way planar data and then scaling the ou-axis value 
to unity. The 2 MHz fundamental, 4 MHz second har- 
monic, and 6 MHz third harmonic average radii are shown 
in Fig. 10a. The 2, 4, and 6 MHz fundamentals are shown 
in Fig. 10b. The average two-way sidelobe levels of the 
second and third harmonic can be seen to be significantly 



lower than those of any of the fundamental harmonics. 
In Fig. 10c, the two 1 MHz average two-way profiles are 
shown. The 4 MHz fundamental curve has a slightly nar- 
rower mainlobe (9.5% at 20 dB down from peak) and sig- 
nificantly higher sidelobe levels than the 4 MHz second 
harmonic profile. Likewise, Fig. 11(d) shows the 6 MHz 
fundamental and 6 MHz third harmonic two-way profiles. 
The 6 MHz fundamental curve is 10.9% narrower at the 
-20 dB level and also shows higher sidelobe levels than 
the third harmonic curve. 

The results shown in Fig. 10 suggest that the second 
and third harmonics maintain lower sidelobe levels than 
the corresponding fundamental harmonics in propagating 
through abdominal wall. In Fig. 11, radial results obtained 
by averaging across five abdominal wall-jittered propaga- 
tions are shown. In Fig. 11a depicts the normalized ra- 
dial averages from the five two-way planar amplitude data 
sets. The average 2 MHz result is overlaid with the corre- 
sponding 4 MHz second harmonic and 4 MHz fundamen- 
tal curves. The -20 dB width of the 4 MHz fundamental 
profile is 6.6%) narrower than the 4 MHz second harmonic 
profile. The -20 dB width for this and subsequent multiple 
propagation-averaged results are given in Table HI. 

Each of the two-way average profiles shown in Fig. 11a 
was then radially summed. Each of the average profiles 
consisted of 499 radial position-magnitude value pairs, 
(n, TTii), i =: 1,. . . ,499. The 499 pairs discretely described 
the average of five abdominal wall-jittered, two-way pro- 
files over a radial extent of 1.2 cm. The discrete radial 
summation of the average profiles was defined by 499 radial 
position-summation value (r^, Sj) pairs, where each sj was 
defined by Sj - tt^Li [(r^-r^.J x (m, + m,_i) /2] . 
The first term in the summation involves the (ro, mo) on- 
axis magnitude value. 

The resulting radial summation or integration profiles 
are shown in Fig. lib. Each of these integrated two-way 
profiles was scaled such that the value at an off-axis radial 
distance of 1.2 cm was unity. The elevated sidelobes of the 
2 and 4 MHz fundamental profiles cause their summation 
profiles to rise significantly beyond the radial extent of the 
mainlobe. This additional rise represents the potential for 
scattering from the sidelobes to significantly reduce the 
contrast resolution of the image. In Table III the radial 
extent at which these integration profiles reach the 0.9 level 
is given. This radial extent is a measure of the sidelobe's 
potential to reduce the contrast resolution of an image. In 
this case the second harmonic's radial extent is 38% less 
than that of its 2 MHz fundamentaL and 63% less than 
the corresponding 4 MHz fundamental radial extent. 

Alternatively, the summation profiles depicted in 
Fig. lib offer the percent of the two-way field's amphtude 
inside or outside a given radius. For example, 91./% of 
the 2 MHz fundamentars amplitude, 96. of the 4 MHz 
second harmonic's amplitude, and 88.6%) of the 4 MHz 
fundamental's amplitude fall inside a radius of 0.25 cm. 
The corresponding percentages falling outside of 0.25 cm 
are 8.3%, 3.3%, and 11.4%, respectively. Ratios of these 
outside percentages could be useful for inferring the rela- 
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live contrasts offered in imaging a low scattering or void 
region of a given size. For a void region approximately 
0.5 cm across, the 4 MHz second harmonic of this device 
might then provide 2.5 times (8 dB) higher contrast than 
the 2 MHz fundamental, and 3.5 times (11 dB) higher con- 
trast than the corresponding 4 MHz fundamental. 

The previous two sets (at 2 and 4 MHz) of five abdom- 
inal wall-jittered propagations were repeated at twice the 
source frequencies. Fig. 12a-b depicts the corresponding 
averaged results from these 4 MHz nonlinear and 8 MHz 
linear propagations. In Fig. 12a the radial average ampli- 
tude curves obtained by averaging the five two-way data 
sets are shown. Note in Fig. 12a that the 8 MHz fundamen- 
tal mainbeam is broader than the corresponding 4 MHz 
fundamental mainbeam. The jitter-imposed lateral resolu- 
tion limits have been encountered, and in fact exceeded, 
at this p6int. Also note that the 8 MHz second harmonic 
mainbeam is narrower than either of the fundamentals'. 
Fig. 12b depicts the additional sidelobe corruption of the 
fundamental beams, and the corresponding increase in the 
second harmonic's relative contrast resolution potential. 

Fig. 13a-b shows average results from five propaga- 
tions through breast wall delay data. In Fig. 13 the re- 
sults considered were from 2 MHz nonlinear and 4 MHz 
linear propagations. In Fig. 13a, the 4 MHz mainbeams are 
significantly broadened over the corresponding abdominal 
wall-jittered mainbeams in Fig. 11a. The sidelobe levels in 
Fig. 13a are also higher than those in Fig. 11a. As was ob- 
served in [7], the breast wall layers produced significantly 
more distortion than the abdominal wall layers. The sec- 
ond harmonic 4 MHz profile is 8.4% narrower at the -20 
dB level than the fundamental 4 MHz. In Fig. 13b the in- 
tegrated profile of the second harmonic has a radial extent 
at the 0.9 level, which is 48% less than the 2 MHz fun- 
damental's and 70% less than the 4 MHz fundamental's. 
Thus, in the more distorting breast wall-jittered propaga- 
tions, the relative advantages of the 4 MHz second har- 
monic were larger than in the abdominal wall propaga- 
tions and included a slightly narrower mainbeam. The 4 
and 8 MHz results for breast wall-jittered propagations 
followed closely the developments previously seen in the 4 
and 8 MHz curves of Fig. 12a-b. Results from these prop- 
agations are included in Table III. 

Finally, in all of the jittered propagations considered, 
the second harmonic mainbeam was narrower than the fun- 
damental mainbeam. The limits on the lateral resolution 
of the linear harmonics though, eventually put limits on 
the second and other higher harmonics. 

VI. Discussion .and Conclusions 

The previous sections showed the predicted liver-path 
beam patterns for a focused Gaussian-apodized transducer 
operating at 2, 4. and 8 MHz. The non- phase aberrated 
propagations suggested that the second and higher har- 
monics formed through finite amplitude distortion should 
have much lower sidelolK; le\'els than their fundamental 



harmonic or the corresponding linear fundamentals. The fi- 
nite amplitude production of these higher harmonic beams 
should allow this sidelobe relationship to hold for any fo- 
cused or unfocused transducer. Pulse propagation analysis 
suggested that the higher harmonics formed in a propa- 
gating pulse-beam can be very well described by consid- 
ering the harmonics produced in the corresponding con- 
tinuous wave propagation. Modeling results also predicted 
that second harmonic levels sufficient for imaging purposes 
could be obtained within the field amplitude limits of the 
mechanical index. 

The introduction of phase jitter as computed from mea- 
sured propagation delays from slices of abdominal wall and 
breast wall caused the sidelobe levels of the second har- 
monic and fundamental beams to rise significantly. In all 
the aberrated propagations considered, the two-way pro- 
file of the second harmonic offered narrower —20 dB main- 
lobe widths and lower sidelobe levels than the fundamental 
beam which produced it. These same second harmonic pro- 
files had slightly broader mainlobes at 4 MHz in abdominal 
wall-jittered propagations than the 4 MHz fundamental 
profiles. In all other jittered propagations considered, the 
second harmonic offered slightly narrower mainlobes than 
the same-frequency fundamental, and substantially lower 
sidelobe levels. Thus, it appears that second harmonic- 
based ultrasonic images could offer significant improve- 
ment in the lateral component of contrast resolution. 

One obstacle to obtaining such images would be arti- 
facts from source contributions to the higher harmonic 
bandwidths. Such artifacts are evident, for example, in 
measured pulses [27]. Perhaps it would not be possible to 
produce source waves with as little energy content outside 
the main bandwidth as that depicted in Fig. 5. 

Fig. 14a depicts such an imperfect source wave. This on- 
axis waveform is depicted in particle velocity units and cor- 
responds to a peak pressure of one half of that of the pulse 
shown in Fig. 5a. The pulse in Fig. 14a heis the same form 
as that in Fig. 5a. with the exception of the initial zero 
portion. Fig. 14b shows the corresponding source pulse 
spectrum. Significant energy content outside the 2 MHz 
bandwidth is visible. The resulting computed focal wave- 
form's spectrum is depicted in Fig. 14c and can be com- 
pared to the earlier focal spectrum in Fig. 5b. Considering 
the log-scaled depiction in Fig. 14c, it is apparent that a 
simple high pass filtering of the received spectrum would 
contain significant contributions from the source. 

Fig. 14(d) shows the resulting nonlinear distortion pulse 
obtained by high-pass filtering the spectrum of Fig. 14c. 
This pulse can be compared to the corresponding distor- 
tion pulse shown in Fig. 5c. In both cases the cut-off fre- 
quency was 3 MHz. In Fig. 14(d) the full 8 microsecond 
period of the computed pulse is shown so that the ringing 
associated with the linear or source content within the sec- 
ond harmonic bandwidth can be seen. This ringing could 
substantially hurt the axial component of the contrast 
resolution and might negate the lateral gains associated 
with the lower sidelobe levels. In this case, improvement 
could be obtained by increasing the source amplitude so as 
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to boost the second harmonic bandwidth levels. Such an 
amphtude-dependent strategy might not be desirable, and 
is further complicated by the fact that off-axis the second 
harmonic drops off much faster than the source content m 
the same bandwidth. 

hi order to eliminate such [)robicms, a two-pulse scheme 
appears to be a good solution. In this method, two source 
pulses would be sent in place of a single pulse in the image 
formation cycle. The two pulses would be identical in form, 



but one would be significantly lower in amplitude. The re- 
ceived echo from this lower amplitude pulse could then 
be used to remove the linear content from the high am- 
plitude pulse. This could be accomplished by subtracting 
an appropriately-scaled version of the received low anii:»h- 
tude signal from the corresponding high amplitude signal. 
The resulting difference signal should then be higli-pass fil- 
tered, followed by the normal sequence of image formation 
steps! The high-p<\ss filtering would be necessary here since 
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piiisc analysis has revcfilod iliat t.lie low frequency content 
of the difference signal is radially wide-spread and conld 
greatly reduce the higiier harmonic sidelobe advantage. 

Fig. 15 shows results relexant to an implementation of 
the two-pulse scheme using the pulse shown in Fig. 1-ia 
a.s the high amplitude source pulse. In Fig. 15a, the focal 
spectrum from 14c is overlaid with the focal spectrum pro- 
(hiced by a half-amplitude version of the same source. The 
low or half-amplitude spectnni) was multiplied by two and 
subtracted from the high amplitude spectrum to obtain the 
difference spectrum depicted in Fig. 15b. Note that if the 
propagations had not involved finite amplitude distortion, 
then this difference spectium would have been all zero. 
The difference spectrum was then high- pass filtered and 
inverse transformed to obtain the effective on-axis distor- 
tion imaging pulse shown in Fig. 15c. Compare this pulse 
with that shown in Fig. 14(d). The transition for the high- 
pass filtering was at 2,75 MHz. This same frequency was 
also appropriate for filtering off-axis difference spectra. 

Initial analysis suggests that this two- pulse scheme ap- 
pears to be capable of extracting the desired largely-second 
harmonic images from any realistic ultrasonic imaging 
pulse. The third harmonic bandwidth depicted in Fig. 15b 
though, does not appear to be separable from the sec- 
ond harmonic bandwidth. The nodal depth between these 
harmonics is not deep enough. Thus, third harmonic (or 
largely-third harmonic) images do not appear to be ob- 
tainable with this two pulse scheme. Unless source pulses 
can be produced with very little energy outside the main 
bandwidth (see the spectrum in Fig. 5b), the production of 
third harmonic images might not be constructive in light 
of axial resolution problems. Initial analysis also suggest 
that a two pulse imaging scheme with a 1/4000 second 
interval between the respective high and low amplitude 
source pulses, should not have significant artifacts due to 
tissue or transducer motion. 

The use of the second harmonic (plus a little contribu- 
tion from higher harmonics) to form an image is an inde- 
pendent alternative to the phase correct ion- based schemes 
which have been and are being examined by other investi- 
gators as a means for improving the contrast performance 
of biomedical ultrasonic imaging. If one or more of these 
phase or beam correction schemes proves successful, then 
they could also be used to improve the contrast of nonlin- 
ear images. To the extent that a phase correction method 
improves the in vivo focusing of a beam it should de- 
crease the beam's sidelobe levels. The corresponding side- 
lobe level decrease for the constituent finite amplitude har- 
monics should be larger than for the fundamental.' 

The higher harmonics may also offer some additional 
opportunities for correcting for beam distortion. The ani- 
phtude of the third harmonic in the focal region is strongly 
affected by the fundamentars amplitude. As discussed in 
Section IV, the nonlinear production of the third harmonic 
is proportional to the product of the amplitudes of the fun- 
damental and second harmonics. This means that much of 
the third harmonic's production occurs in the focal region 
where beam distortion can significantly reduce the am- 
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Pig. 15. The log-scaled focal spectrum from 14c overlaid with the 
focal spectrum from the same source using a half amplitude versioti 
of the source pulse (as depicted in Fig. 14a). (b) The resulting differ- 
ence spectrum computed for the two spectra shown in (a), (c) The 
corresponding nonlinear distortion pulse obtained by constructing 
with the spectral information in (b) starting at 2.75 MHz. 



plitude of the fundamental and second harmonic. It thus 
might be possible to develop an iterative scheme to correct 
for beam de foe using, using the amplitude of the received 
tliird harmonic for feedback. 

Second harmonic images may also provide for a means 
of reducing speckle. In particular, an image formed as the 
sum of a second liarmonic image and the corresponding 
fundamental image should have less speckle than either of 
the constituent images. Since the second harmonic is twice 
the frequency of the fundamental and has a largely con- 
stant phase relationship with the fundamental in the main- 
lobe, the second harmonic image's speckle pattern should 
be conveniently out of phase with that of the fundamental. 



f 



CHRISTOPIIKR: KINITE AMPLtTi;iJI:: OISTOTTTION-BASKO I N I lOMOG KN EOUS PULSK ECHO ULTRASONIC IMAGING 139 



Finally, statistically-defined jitter planes were also used 
(as an alternative to phase delay planes directly computed 
from measured temporal delays). Using the RMS temporal 
fluctuation and spatial correlation lengths computed for 
the abdominal wall [4] and breast wall [7] measurements 
and the algorithms given in (28], jitter planes were com- 
puted. Unfortunately, these simulated jitter planes gave 
markedly different defocusing results than the meixsured 
delay planes, so they were abandoned. 

The results of this study suggest that the finite am- 
plitude distortion of biomedical ultrasonic imaging pulses 
could be utilized to enhance the contrast resolution of im- 
ages. Experimental work is planned- Additionally, other 
uses for finite amplitude distortion may be possible as well. 

Acknowledgments 

Helpful discussions with E. L. Carstensen, X. Chen, L. 
M. Hinkleman, K. N. Kumar, D. L. Liu, T. D. Mast, K. J. 
Parker, and R. C. Waag are gratefully acknowledged. The 
efforts of L. M. Hinkleman and D. L. Liu to provide the 
measured tissue delay data were invaluable. 



References 

[1] M. Moshfeghi and R. C. Waag, "In vivo and in vitro ultrasound 
beam distortion measurements of a large aperture and a con- 
ventional aperture focused transducer," Ultrasound in Med. and 
Biol, vol. 14, pp. 415-428, 1988. 

[2) L. Nock and G. E. Trahey, "Phase aberration correction in med- 
ical ultrasound using speckle brightness as a quality factor," J. 
Acousi. Sac. Am., vol. 86, pp. 1819-1833, 1989. 

[3j S.H.P. Bly, D. Lee-Chahal, D. R. Foster, M. S. Patterson, F. S. 
Foster, and J. W. Hunt. "Quantitative contrast measurements 
in S-mode images comparison between experiment and theory," 
Ultrasound in Med. and Biol., vol. 12, pp. 197-208, 1986. 

(4) L. M. Hinkleman, D. L. Liu, L. A. Metlay, and R. C. Waag, 
"Measurements of ultrasonic pulse arrival time and energy level 
variations produced by propagation through abdominal wall," 
J. AcousL Soc. Am., vol. 95. pp. 530-541, 1994. 

[5] D. L. Liu and R. C. Waag, 'Time-shift compensation of ultra- 
sonic pulse focus degradation using least-mean-square estimates 
of arrival time," J. Acoust. Soc. Am., vol. 95, pp. 542-555, 1994. 

[6] D. L. Liu and R. C Waag, 'Correction of ultrasonic wave- 
front distortion using back propagation and a reference wave- 
form method for time-shift compensation," J. Acoust. Soc. Am., 
' vol. 96, pp. 649-660, 1994. 

(7) L; M. Hinkleman, D. L. Liu. L. A. Metlay, and R. C. Waag. 
"Measurements of ultrasonic pulse arrival time and energy level 
variations produced by propagation through breast wall," J. 
Acoust. Soc. Am., vol 95, pp. 530-541, 1994. 

(8) R. Mallart and M. Fink, ".Adaptive focusing in scattering media 
through sound speed inhomogeneities: The van Citlert Zernickc 
approach and focusing criterion." J. Acoust. Soc. Am., vol. 96, 
pp. 3721-3732, 1994. 

(9| Q. Zhu and B. D. Steinberg, ' VVavefront amplitude distribution 
in the female breast," J. Acoust. Soc. Am., vol. 96, pp. 1-9, 
1994. 

[10] Q. Zhu and B. 0. Steinberg, -Wavefront amplitude distortion 
and image sidelobe levels: Part 1 - Theory and computer simula- 
tions." IEEE Traris. Ullrason. FeiToelec. Freq. Contr., vol. 40, 
pp. 7^17-753, 1993. 

[11] Q. Zhu and B. D. Steinberg, "Wavefront amplitude distortion 
and image sidelobe levels: I^ut W - In vivo experiments," IEEE 
Trans. Ulirasoii. Ferroelec. Freq. Conir.,. vol. 40, p|). 747-753, 
1993. 



[12] J. Y. Lu, H. Zou, and .L F. Grecnleaf, "Biomedical ultrasound 
beam forming," Ultrasound in Med. and Biol., vol. 20, pp. 403- 
428, 1994. 

[131 Karaman, A. Atalar, H. Koymen, and M. O'Donnell, "A 
phase aberration correction method for ultrasound imaging," 
IEEE Trans. Ullrason. Ferroelec. Freq. Contr., vol. 40, pp. 275- 
288, 1993. 

[14] .J. A. TenCate, "An experimental investigation of the nonlinear 
pressure field produced by a plane circular piston," J. Acoust. 
Soc. Am., vol. 94, pp. 1084-1089, 1993. 

[15) A. C. Baker, K. Anastasiadis, and V. F. Humphrey, "The non- 
linear pressure field of a plane circular piston; theory and ex- 
periment," J. AcousL Soc. A7n., vol. 78, pp. 1483-1487. 1988. 

(161 'T- Christopher and K.J. Parker, "New approaches to nonlin- 
ear diffractive field propagation," J. Acoust. Soc. Am., vol. 90, 
pp. 488-499, 1991. 

(17) M. A. Averkiou and M. F. Hamilton, "Measurements of har- 
monic generation in a focused finite-amplitude sound beam," J. 
Acoxist. Soc. Am., vol. 98, pp. 3439-3442, 1995. 

[18) T. Christopher, "Modeling acoustic field propagation for medi- 
cal devices," Ph.D. thesis, U. of Rochester, 1993. 

[19] T. Kamakura, M. Tani, Y. Kumamoto, and K. Ueda, "Harmonic 
generation in finite amplitude sound beams from a rectangular 
aperture source," J. Acoust. Soc. Am., vol. 91, pp. 3144-3151, 
1992. 

[20] A. C. Baker, A. M. Berg, A.- Sahin, J. N. Tjotta, "The nonlinear 
pressure field of plane, rectangular apertures: Experimental and 
theoretical results," J. Acoust. Soc. Am., vol. 97, pp. 3510-3517, 
1995. 

[21] M. E. Lyons and K. J. Parker, "Absorption and attenuation in 
soft tissues II - experimental results," IEEE Trans. UUrason. 
Ferroelec. Freq. Contr., vol. 35, pp. 511-521, 1988, 

(22) S. A. Goss, R. L. Johnston, and F. Dunn, "Compilation of 
empirical ultrasonic properties of mammalian tissues. II.," J. 
Acoust. Soc. Am., vol. 68, pp. 93-108, 1980. 

[23] W. K. Law, L. A. Frizzell. and F. Dunn, "Determination of the 
nonlinearity parameter B/A of biological media," Ultrasound in 
Med. and Biol., vol. 11, pp. 307-318, 1985. 

[24] P. T. Christopher, "A nonlinear plane wave algorithm for 
diffractive propagation involving shock waves," J. Comp. 
AcousL, vol. 1, pp. 371-393, 1993. 

[25] T. Christopher, "Modeling the Dornier HM3 lithotripter," J. 
AcousL Soc. Am., pp. 3088-3095, 1994. 

[26] J. W. Goodman, Introduction to Fourier Optics. New York: 
McGraw-Hill, 1968. 

[27] A. C. Baker and V. F. Humphrey, "Distortion and high- 
frequency generation due to nonlinear propagation of short ul- 
trasonic pulses from a plane circular piston," J. AcousL Soc. 
Am., pp. 1699-1705, 1992. 

[28] D. L. Knepp, "Multiple phcise-screen calculation of the temporal 
behavior of stochastic waves," Proceedings of the IEEE 71(6), 
pp. 722-737, 1983. 



Ted Christopher was born in Baltimore, Maryland in 1960. He 
received the B.S. degree in Computer and Information Science from 
the University of Massachusetts, Amherst, MA in 1981 and the M.S. 
and Ph.D. degrees in Electrical Engineering from the University 
of Rochester, Rochester. NY in 1984 and 1993, respectively. He is 
currently working at the Diagnostic Ultrasound Research Labora- 
tory at the University of Rochester. He has held a number of ul- 
trasound research and teaching- related positions at the University 
of Rochester over the last 12 years. His research interests include 
modelling biomedical ultra.sonic wave propagation and improving the 
quality of biomedical nltriLSonic images. 



2. 



t 

r«M...iUrv ned^>^^nn nf Sharon MuKagh, M B^ O | rn 

I, SharonL.Mulvagh,M.D., declare thatlbeUeve the following facts to be o^^^^ 
peisonal knowledge: CD 
1 . I am a medical doctor on staff a. the Mayo M«Ucal Clime, and have held this M<ioB since 
November 1989. 

Expetiroental and CUnical Studies, New Technologies. 

3 My presentation inclnded a videotape shov«ng the ^ff '""^ ?J 

video images. 

4 -n^ videotape I presented showed images of adog's beattag heart, in ^^^^"'^ 

followed by hamomc images ofthe contrast agent, incs/^i'T-. 
represents the datt on which the taped expethnent was conducted. 

I have enclosed two videotapes. The ^V^^^}^^^^ "^iS^d vid^tape is 

9-28-94, Complete Copy". 

6. The confaenc. presentadon was open to all who '^•^'^^'^^L. to be 
confidentiaUty on the part of the attendees. Accordingly, I considered my oisci 

public. 

, h«,by declare under penalty of penury that the foregoing U, to the best of my knowledge, 

correct. 

Dated this of MmO^ ' 1998, by 
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Evidentiary Declaration of Gregory L. Hollev S 

O ^ O 
^ ^ 

I, Gregory L. HoUey, declare that I believe the following facts to be correct based on^jOiy^rsdn^J 
knowledge: 

^ ^ o 

1 . I am an Image Analysis Engineer at Acuson Corporation, and have held this posi^S^n since 
July 20, 1987. I have been involved in the development of contrast imaging and second 
harmonic imaging at Acuson for the last several years. 

2. I have reviewed the videotape referred to in the Declaration of Sharon L. Mulvagh, M.D. 

3. The videotape shows that the Acuson 128 system was transmitting ultrasound energy into a 
dog's beating heart at a fundamental frequency and receiving ultrasound echoes at a harmonic 
frequency, irrespective of whether or not contrast agent had been injected into the dog. 

4. The tape shows harmonic imaging, in the absence of contrast agent, in the section entitled 
"Baseline Imagent". During this part of the imaging, a second harmonic image of the tissue 
of the dog's heart is shovm. For approximately the first eight seconds of this image, no 
contrast agent is visible in the image. The image that is visible is formed at least in part from 
the second harmonic response of the tissue to the transmitted ultrasound wave. The image is 
not formed from harmonic scattering by a contrast agent. 

5. I have reviewed the paper by Ted Christopher entitled, "Finite Amplitude Distortion-Based 
Inhomogeneous Pulse Echo Ultrasonic Imaging," (IEEE Trans. USFC 44(1) (Jan. 97) pp. 
125-139). Finite amplitude distortion during the propagation of the ultrasound wave through 
the body, as described in Christopher's paper, is the mechanism by which the second 
harmonic response of the tissue is generated. This response contributes to the displayed 
image of the dog's heart that is based upon echoes received at the harmonic frequency. 

6. I also reviewed the videotape entitled "Contrast Echocardiography Investigations," referred 
to in paragraph 13 of Joan C. Main's declaration. 

7. In a segment of this videotape entitled "Experiment 1", a video clip is shown on the right side 
of the screen labeled "Second harmonic with contrast." This video clip shows a short axis 
view of an animal's heart. The first second or so of this clip shows a second harmonic image 
of the animal's heart in the absence of contrast agent. The display shows a second harmonic 
image of the tissues of the heart itself. This image derives at least in part from second 
harmonic distortion that occurs as the transmitted ultrasound pulse propagates through the 
tissues of the animal subject. In my opinion, such second harmonic distortion makes a 
measurable contribution to these images. As noted in paragraph 5, "finite amplitude 
distortion" is a name for the mechanism by which this second harmonic distortion is 
generated. 



8. 



Another example of images formed in part from second harmonic generated by finite 
amplitude distortion is shown in the segment labeled "Experiment 2." Again, the right side 
of the screen shows a second harmonic video clip. For at least the first several frames, there 



is no contrast agent visible in the left ventricle of the heart (the image is oriented so the left 
ventricle is visible on the right side of the video clip). In these first several firames, the image 
of the tissues of the left ventricle seen in the clip is formed at least in part from second 
harmonic distortion of the transmitted pulse caused by propagation through the tissues of the 
animal, and not from second harmonic scattering from the contrast agent. In my opinion, 
such second harmonic distortion makes a measurable contribution to the images in these first 
several frames. 



I hereby declare under penalty of perjury that the foregoing is true and correct 



Dated this fj_ day of MoAx^i^ ^ 1993, by 
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Evidentiary Declaration of Joan C. Main 

I, Joan C. Main, declare that I believe the following facts to be correct based on my 
personal knowledge: 

1 . I have been employed in the marketing department of Acuson Corporation s^e Mgy ^ 
1986, My primary focus over the last few years has been in the area of contr^ ^ ^ 
imaging. In connection with my responsibilities in this area, I have worked claSely g 
with physicians having a research interest in the field. up 

Cm 

2. I also worked with an Acuson engineer, Paul Chandler, in supporting second hafiaionic 
engineering development efforts. He made modifications to our Acuson 128 
ultrasound system to support transmission at a fijndamental center frequency of 
2.5MHz and reception at a harmonic center frequency of 5.0MHz. Special wideband 
transducers were developed for these investigations by Acuson' s Principal Fellow, 
Amin Hanafy. 

3. Among the initial physicians with whom I worked were Doctors Sharon Mulvagh and 
Tony DeMaria. In working with these physicians, I assisted in the operation of the 
modified Acuson 128 system, which is internally known and labeled as System 1 122. 

4. Specifically, in working with Doctor Sharon Mulvagh in August 1994, 1 brought 
Acuson' s System 1 122 to the Mayo Clinic to do dog studies with the Imagent® 
contrast agent from Alliance Corporation in San Diego. We videotaped the dog's 
beating heart while imaging with System 1 122, in the harmonic mode. Baseline 
images were gathered first. Then contrast agent was injected into the dog's vein, and 
contrast images were gathered. I personally did at least some of the imaging for this 
study. 

5. Segments from this study were shown in a presentation that Dr. Mulvagh made and 
that I attended, on or about September 30, 1994, in the Knickerbocker Hotel in 
Chicago, Illinois, at an echocardiography conference entitled "Advances in 
Echocardiography: Contrast Echocardiography Perfijsion Imaging Transesophageal 
Echo." Her presentation was part of Session I of the day. The Session was entitled, 
"New Agents: Experimental and Clinical Studies, New Technologies." 

6. Her presentation included a videotape showing the use of the contrast agent Imagent® 
on an Acuson 128 machine. Images on this tape came from the study that we worked 
on together. The tape is identified by Mayo Clinic ID#CSI, and Echo Date 9-28-94. 
The machine was configured to transmit ultrasound energy at a fijndamental frequency 
(2.5MHz), and receive ultrasound echoes from the anatomy and/or contrast agent at a 
harmonic (5MHz) of that frequency. This configuration is indicated by the highlighted 
5.0MHz icon on the video images. 

7. The videotape that she presented showed images of a dog's beating heart, taken on 
8/22/94, both in the absence of contrast agent (called "Baseline Imagent" on the tape), 
as well as with contrast agent subsequently injected into the dog. In her presentation, 
she pointed out the baseline image prior to the contrast injection. Accordingly, 




harmonic images of the tissue were seen, followed by harmonic images of the contrast 
agent. 

8. The conference presentation was open to all who attended, with no obligations of 
confidentiality on the part of the attendees. 

9. Turning now to my earlier work with Dr. DeMaria, in early 1993 I accompanied Paul 
Chandler and his clinical marketing team to conduct a proof of principle experiment. 
We did a pig and a dog study using System 1 122 at the University of California at San 
Diego, with the Schering contrast agent called SSH508 (now known as Levovist®). 

10. The studies were documented on Betacam videotape, and segments were chosen and 
edited for a private showing to Acuson's key luminaries worldwide. Paul Chandler 
and I presented the information on the tape to this group of physicians under a non- 
disclosure agreement on Sunday, March 14, 1993, prior to the start of the American 
College of Cardiology meeting in Anaheim, California. 

1 1. On June 7, 1993, 1 sent the tape to TVA Incorporated in Mountain View to be 
converted to Super- VHS format for additional presentations to physicians and their 



12. Later in June of 1993, 1 personally presented the Betacam version of the videotape to 
key investigators on the show floor at the American Society of Echocardiography. 
This was a public showing, in that no non-disclosure agreements were signed or 
requested. 

13. Enclosed is the Super- VHS version of the same videotape. It has the same images as 
were on the Betacam version that I presented at the American Society of 
Echocardiography. 

14. This tape showed two experiments. In experiment number one, on the right side of the 
screen is a short axis view of an animal heart. For the first second or so of this video 
clip, there is no contrast agent present; the display just shows a harmonic image of the 
tissues of the heart itself 

15. The second experiment shows an apical 4-chamber view showing a second harmonic 
image of the heart while contrast agent is injected. For the first few fi-ames, there is no 
contrast agent present in the left ventricle of the heart (on the right side of the image). 
The image of the left ventricle of the heart is formed at least in part from the second 
harmonic response of the tissue itself and not from a contrast agent. 

I hereby declare under penalty of perjury that the foregoing is, to the best of my 
knowledge, correct. 



staff. 



Dated this 1 1 day of February 1998, by 




